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DUAL ACCELEROMETER USAGE STRATEGY FOR ONBOARD 

SPACE NAVIGATION 

Renato Zanetti and Chris D’Souza f 


This work introduces a dual accelerometer usage strategy for onboard 
space navigation. In the proposed algorithm the accelerometer is used 
to propagate the state when its value exceeds a threshold and it is used 
to estimate its errors otherwise. Numerical examples and comparison to 
other accelerometer usage schemes are presented to validate the proposed 
approach. 

INTRODUCTION 

In the presence of significant and consistent non-gravitational accelerations, accelerom- 
eter measurements are often used in lieu of analytical expressions to propagate the state of 
aerospace vehicles in model-based estimation algorithms such as the Kalman filter [1, 2]. 
In space applications, accelerometers are often used to propagate the state of the vehicle 
through a translational maneuver, or “bum”. At all other times non-gravitational forces are 
usually much smaller than the accuracy of most commercially available accelerometers, 
therefore space vehicles commonly threshold the acceleromter, i.e. they use it only when 
the measurements are above a predetermined value [3]. The Space Shuttle rendezvous and 
proximity operations program (RPOP), for example, employs this strategy [4]. 

Rarely accelerometers are used as external measurements to update the state of the ve- 
hicle. One example of this approach is Ref. [5] in which accelerometers are used in con- 
junction with a filter bank for Mars entry navigation. Under this type of implementation, a 
model of the non-gravitational forces is required. When on a launch pad a vehicle is sta- 
tionary with respect to Earth and the predicted accelerometer measurement is very easily 
obtained. Pre-launch operations therefore commonly employ accelerometer measurements 
to update the estimate of its repeatable errors, such as biases, misalignments, and scale 
factors. 

Accelerometers usually have two types of repeatable biases. A constant bias (sometimes 
represented as a first order Markov process [6] with very large time constant of the or- 
der of 24 hours) and a faster changing Markov process with time constant of around one 
hour. Velocity random walk (VRW) also corrupts the measurement, but this source of error 
is not estimable because it is white acceleration noise. In the standard accelerometer us- 
age by space vehicles the accelerometer errors are estimated from external measurements 
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through their correlation to other states. This correlation is built during maneuvers at which 
times the accelerometer measurement is incorporated into the filter dynamics. In between 
maneuvers the accelerometer measurement is not used at all, therefore the estimate of the 
one-hour Markov bias degrades, which results in a worst navigation performance during the 
next maneuver. For vehicles with small thrusters and maneuvers far apart, this approach 
is often not sufficient to meet mission goals. A common practice for these vehicles is to 
use a simple averaging scheme to estimate the accelerometer bias prior to maneuvers. This 
average is carried outside the navigation filter, therefore no correlation between the filter 
error and the accelerometer bias estimation error is taken into account when formulating a 
navigation solution. 

This work introduces a dual accelerometer usage strategy for onboard space navigation. 
In the proposed algorithm the accelerometer is used to propagate the state when its value 
exceeds a threshold and it is used to estimate its errors otherwise. The paper is organized 
as follows: accelerometer thresholding is first introduced, the navigation algorithm is then 
presented and validated with numerical simulation, finally conclusions are drawn. 

ACCELEROMETER THRESHOLDING 

This section presents a simple example to show the need of accelerometer thresholding 
in spacecraft navigation. Outside of thrust, drag is usually the biggest non-gravitational 
acceleration source in low Earth orbit (LEO). The international space station (ISS) is used 
as an example to quantify the contribution of drag. The ISS is a large structure placed 
at an altitude of 400 km and is subject to high drag. A 2008 analysis from the Euro- 
pean Space Agency* finds that the ISS coefficient of drag is 2.07, the frontal area between 
700 and 2300m 2 depending on the configuration, and the yearly average air density 3.98 
10 _12 kg/m 3 . Assuming a configuration with frontal area 1000m 2 the drag acting on the 
ISS is 0.25N. The mass of ISS was 2.5 10 5 kg in 2008, resulting in a drag acceleration of 
10~ 6 m/s 2 or OA/ig. 

A relatively good accelerometer can have a one hour markov with steady-state standard 
deviation a ss = 10 fig. For each accelerometer’s axis, the markov bias b a evolves as 

1 

b a = --b a + v , (1) 

r 

where r = 3600s and v is a zero-mean white process with spectral density S v = o 2 ss j (2 r). 
In order to establish to which accuracy this markov bias can be estimated, we assume a 
continuous Kalman filter, we also assume no non-gravitational forces acting on the vehicle, 
therefore the accelerometer measurement is 

y = b a + n (2) 

for each axis. The velocity random walk rj is a zero-mean white process, a spectral density 
S v = (10 fig'/s) 2 is assumed, which is a relatively good value. The estimate of the markov 

*http : / / www . esa . int / esaMI/ Space_In_Bytes/ SEM7FXJ2 6DF_0 . html 
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bias evolves as 


( 3 ) 


A 1. „ 

ba = b a + K(y - b a ), 

T 

where K = P/S v and P is the estimation error variance evolving as 

P = ~-P + -<t 2 ss -P 2 /S v . (4) 

r r 

After some algebra it follows that the steady-state value of the estimation error standard 
deviation is \/P ss = 1.5//g. Even under the optimistic assumptions of this example the 
estimation error of the markov bias is 15 times bigger than drag. When using the com- 
pensated accelerometer measurement to propagate the state an additional error source is 
the VRW. During non-thrusting phases it is therefore much more accurate not to use the 
accelerometer to propagate the state but to use a simple drag model. 


DUAL ACCELEROMETER USAGE STRATEGY 


In all spacecraft missions the authors are aware of accelerometers are thresholded during 
orbital coast flight (sometimes this operation is referred to as accelerometer gating). In 
the majority of this mission the accelerometer bias is estimated to improve the navigation 
solution. The estimation occurs in either of two ways. A common solution is to include the 
accelerometer bias as a state in the filter and to estimate it through external measurements 
and the correlations built during maneuvers. The issue with this approach is that when 
maneuvers occur far apart the estimate degrades in-between maneuvers, often resulting in 
a poor knowledge of the bias when the following maneuver occurs. The usual solution to 
this problem is to have an external estimator of the accelerometer bias. While this sec- 
ond solution usually produces good results, it completely ignores the inevitable correlation 
between the states in the filter and the accelerometer bias estimate. To optimally account 
for this correlation we propose an integrated filter that uses the estimate of the bias during 
maneuvers and estimates it during coast flight. 

Accelerometer measurements fed to the navigation filter are usually integrated accel- 
erations over the last time step and they are compensated for sculling errors. Sculling 
compensation means that the effects of the rotation of the vehicle during the time step are 
compensated and the measurement Av| is an inertial change in velocity between times 
tk- 1 and t k coordinatized in the body-fixed frame at time t k . When the filter is called it first 
propagates the position (r), velocity (v), and accelerometer bias (b) to the current IMU 
time tk as 


d 

d 
d 

dl 

±b i ’ = -lbi 

dt r a 


— r = v 

dt 


** = * + k% 


(5) 

(6) 
(7) 


When the accelerometer measurement is above a threshold the estimated non-gravitational 
acceleration a* is obtained from the accelerometer measurement a* = T l(t k ) (Av^/A t — 
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b 6 ). When the measurement is below the threshold a* is either zero or a model of drag, 
depending on the application. Matrix T \(t k ) is the coordinate transformation from the 
body-fixed frame at time t k to the inertial frame. Other states such as attitude can easily be 
incorporated in the filter. 

When the accelerometer is included during propagation the state is only updated with 
external measurements. When the accelerometer is thresholded an estimate of the mea- 
surement Av b k is formed. This estimate is given by 

r t k nt k 

A\ b k = / (bf t + d b )dt = r(e A */ r — l)b(f fc ) + / d b dt ~ At (b(t k ) + d 5 ), (8) 

''tk - 1 •'tk - 1 

The measurement mapping matrix for the accelerometer measurement, H„ is a 3 x n matrix, 
where n is the number of states in the filter. Matrix H a has zeros everywhere except for an 
identity matrix times At at the 3x3 block corresponding the accelerometer bias state. The 
Kalman gain is calculated as usual 

K t = P^Hj(H,.P,-H T + R„)- 1 , (9) 

where P k is the a priori estimation error covariance matrix and R a is the accelerometer 
noise covariance. The accelerometer white noise is expressed in terms of a velocity random 
walk with an associated spectral density S a whose units are the square of m/s/ a/s or, more 
frequently, pgy/s. Therefore the covariance of the noise over an IMU step At is given by 

Ra = s a At. (10) 

It is not desirable to update all the states, only the accelerometer bias should be updated. 
Therefore a consider gain K( is formed by zeroing all the rows of K/ corresponding to the 
other states. Since the new gain is not optimal the Joseph’s formula is used to obtain the 
a posteriori covariance 

Pi = (i»*„ - K;ii„iP- (i„ xn - kjh„) t + k;.r„(k;) t od 

The state is updated as 

Xfc = Xfc + K£ (Av£ - Av£) . (12) 

The choice of the threshold value is usually driven by the thruster size, the accelerometer 
accuracy, and engineering judgement. A good rule of thumbs is to threshold the accelerom- 
eter when the measurement is below the 3 a value of the sum of the bias and noise, assumed 
independent 

(Av b k ) T Av b k < 9 (trace B a Af 2 + trace S a A t) , (13) 

where B a is the accelerometer bias covariance from the IMU specifications. If Av| is 
compensated with the estimate of the bias, B a is replaced with the accuracy of the bias 
estimate. It is also possible to choose a non-constant threshold 

(Avj( — b 6 At) T (Av^, — h b At) < 9 (trace P m A t 2 + trace S a Af) , (14) 
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where P IM is the 3x3 portion of the covariance corresponding to the accelerometer bias 
state. Spacecrafts usually operate at low angular rates, but in case of high rates the term 
b b At in Eq. (14) needs to be replaced by 



rpfr(tfc) 

A &(t) 


b b 


(r) dr, 


(15) 


where is the transformation matrix that takes the vehicle’s body frame at time r into 
the body frame at time t/-. 

This section assumes that the entire accelerometer measurement is either applied or 
thresholded, it is also possible to develop similar algorithms by considering the measure- 
ment from each axis independently. 

NUMERICAL RESULTS 

An orbital rendezvous is used as an example. The target spacecraft is placed on a circular 
400km orbit. The chaser vehicle is placed in a circular orbit 4km below the target and 
14km behind. At a downrange of 13.2km an altitude rase maneuver is commanded. This 
maneuver is a 135 degree lambert targeting taking the chaser to 1 .4km below and 4.3km 
behind the target. A clean up maneuver is performed half the way during this transfer. 
A circularization maneuver is performed once the end of the transfer (2030s). A third 
maneuver is performed that will take the chaser 500m below the target with zero relative 
velocity. Figure 1 shows the relative trajectory. The origin is the location of the target, 
while the line represent the relative position in time. The y - axis is the altitude, while the 
x-axis is downrange. 


In-Plane LVLH Trajectory 



Figure 1. Relative In-Plane LVLH Trajectory 
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The simulated true environment contains a 9 by 9 gravity model and atmospheric drag. 
The filter is a nine state filter with relative position, velocity, and accelerometer bias. The 
propagation includes J2 effects and no drag. The chaser vehicle has a 20N thruster and a 
mass of 1000kg. The accelerometer errors are given in Table 1. 


Error Type 

1 a Error 

Units 

Time Constant 

3600 

s 

Markov Bias 

100 

Fg 

Velocity Random Walk 

100 

Fgv^ 


Table 1. Accelerometer Model Error Parameters 


Together with the accelerometer there is a long range radar whose errors are shown in 
Table 2. The radar provides range to the target as well as azimuth and elevation to it. 


Error Type 

la Error 

Units 

Range Error 

5 

m 

Angles Errors 

0.5 

deg 


Table 2. Radar Model Error Parameters 


The initial filter covariance is diagonal with entries as shown in Table 3, while the veloc- 
ity process noise covariance is 10 -6 m 2 /s 3 . 


Error Type 

la Error 

Units 

Relative Position 

10 

m 

Relative Velocity 

0.1 

m/s 

Accelerometer Bias 

100 

Fg 


Table 3. Initial Estimation Errors 


Figs. 2 and 3 show the filter’s performance without thresholding the accelerometer (i.e. 
the accelerometer measurement is always used to propagate the vehicle state). The lines 
show the estimation error standard deviation. Figs. 4 and 5 show the filter’s performance 
with thresholding the accelerometer. Figs. 6 and 7 show the filter’s performance when the 
accelerometer bias is estimated from the accelerometer measurements outside maneuvers. 
Figs. 8 and 9 show a direct comparison of the three methods. It can be seen that the new 
algorithm outperforms the simple thresholding scheme during firing of the thrusters. 

CONCLUSIONS 

This papers presents a dual accelerometer usage in an orbital Kalman filter. The ac- 
celerometer is both used to propagate position and velocity during maneuvers and to update 
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RSS Velocity (m/s) 


Spacecraft Position Estimation Error 



Figure 2. Position Estimation Error Without Thresholding 


RSS Spacecraft Velocity Estimation Error 



Figure 3. Velocity Estimation Error Without Thresholding 
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RSS Velocity (m/s) 


Spacecraft Position Estimation Error 



Figure 4. Position Estimation Error With Thresholding 


RSS Spacecraft Velocity Estimation Error 



Figure 5. Velocity Estimation Error With Thresholding 
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RSS Velocity (m/s) 


Spacecraft Position Estimation Error 



Figure 6. Position Estimation Error New Algorithm 


RSS Spacecraft Velocity Estimation Error 



Figure 7. Velocity Estimation Error New Algorithm 
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RSS Spacecraft Position Estimation Error 



Figure 8 Comparison of Position Estimation Error, green line shows no threshold, 
black line with threshold and red line is the new algorithm 


RSS Spacecraft Velocity Estimation Error 



Figure 9 Comparison of Velocity Estimation Error, green line shows no threshold, 
black line with threshold and red line is the new algorithm 
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the state outside of maneuvers. The advantage of this approach is its superior performance 
than simply thresholding the accelerometer. The correlation between accelerometer bias 
and position and velocity during the maneuvers is not sufficient to adequately estimate the 
bias during coast flight. Therefore the estimate of the bias degrades outside of maneuvers 
adding considerable uncertainty during subsequent maneuvers. A common solution to this 
problem is to estimate the accelerometer bias outside of the navigation filter just prior to 
a maneuver is performed. A simple averaging scheme is often used, but a Kalman filter is 
also a possibility. The advantage of the proposed scheme over a separate bias estimator is 
that a single estimator is globally optimal because it accounts for all the correlations. 
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